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An unusual set of orbits about extreme Kerr black holes resides at the Boyer-Lindquist radius 
r = M, the coordinate of the hole's event horizon. These "horizon skimming" orbits have the 
property that their angular momentum L z increases with inclination angle, opposite to the familiar 
behavior one encounters at larger radius. In this paper, I show that this behavior is characteristic of 
a larger family of orbits, the "nearly horizon skimming" (NHS) orbits. NHS orbits exist in the very 
strong field of any black hole with spin a > 0.952412M. Their unusual behavior is due to the locking 
of particle motion near the event horizon to the hole's spin, and is therefore a signature of the Kerr 
metric's extreme strong field. An observational hallmark of NHS orbits is that a small body spiraling 
into a Kerr black hole due to gravitational-wave emission will be driven into orbits of progressively 
smaller inclination angle, toward the equator. This is in contrast to the "normal" behavior. For 
circular orbits, the change in inclination is very small, and unlikely to be of observational importance. 
I argue that the change in inclination may be considerably larger when one considers the evolution of 
' inclined eccentric orbits. If this proves correct, then the gravitational waves produced by evolution 

through the NHS regime may constitute a very interesting and important probe of the strong-field 
nature of rotating black holes. 
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I. INTRODUCTION 



> 
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The space-based gravitational-wave detector LISA |jj is being designed to make very precise measurements of the 
characteristics of black hole spacetimes. One source that is particularly well-suited for such measurements is the 
inspiral of a small (mass \x = 1 — 10 M©) compact body into a massive (M = 1O 5 ~ 7 M0) Kerr black hole. Depending 



upon the values of M, /i/M, and the hole's spin a, such an inspiral will spend several months to several years in 
LISA's frequency band radiating 10 5 — 10 6 gravitational-wave cycles. By accurately measuring these cycles, LISA 
should be able to build a "map" HQ of the spacetime, testing in detail the predictions of general relativity. 

As the community begins developing strategies for analyzing LISA's datastream, it is important to revisit and 
^ ri i carefully analyze the sources one expects to measure. For extreme mass ratio inspirals, this means understanding 
more deeply the character of Kerr black hole orbits and the nature of gravitational radiation reaction deep in the Kerr 
metric's strong field. Of particular interest are features that might simplify data analysis (which is likely to be very 
difficult, given the many cycles that must be tracked) or that might constitute a strong signature of the spacetime. 
In this paper, I analyze what might be such a feature — a unique signature of the inspiral of a body through the 
C$ ' extreme strong field of rapidly rotating black holes. 

The key piece of this analysis was first discussed by Wilkins. Reference || contains a detailed examination of 
circular orbits of extreme (a = M)F] Kerr black holes. (In this context, "circular orbit" means "orbit of constant 
Boyer-Lindquist coordinate radius" .) One very interesting result given in |5| is the existence of "horizon skimming" 
orbits. These are circular orbits of varying inclination angle at the same coordinate radius as the event horizon, 
r = M. (Despite being at the same radial coordinate, one can show that these orbits have distinct proper separation, 
and in particular lie outside the event horizon; see Ref. [^J, particularly Fig. 2.) An extremely interesting feature of 
the horizon skimming orbits is that as an orbit's inclination angle l is increased, its angular momentum component 
L z likewise increases: dL z jdi > for the horizon skimming orbits. This property holds over a sizable range of radius, 
out to r ~ 1.8M. I will call the full set of orbits for which dL z /di > the "nearly horizon skimming" (NHS) orbits. 
This defining property of NHS orbits is opposite to weak-field intuition. For example, in Newtonian theory orbits at 
constant radius have L z = \L\ cos/, (where \L\ is the same for all orbits at radius r), which decreases as t increases. 

Intuition from Newtonian theory is highly suspect in the strong field of black holes. However, the true behavior of 
L z as a function of l is qualitatively the same as in Newtonian theory over a wide range of orbital radii and spins. For 



I use units where G = 1 = c. 
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example, when a = 0.95M, dL z /di < for all circular orbits, so there are no NHS orbits when a — 0.95A/. It turns 
out that NHS orbits can only exist when a > 0.952412Af . This is the smallest spin for which stable orbits come close 
enough to the hole's event horizon that dL z /di can switch sign: the property dL z /di > arises because, very close to 
the black hole, all physical processes become "locked" or "frozen" to the hole's event horizon (7) . In particular, their 
orbital motion locks to the horizon's spin. This locking dominates the "Keplerian" tendency of an orbit to move more 
quickly at smaller radius (fKepier = y/M/r), forcing a body to actually slow in the innermost orbits. The locking is 
particularly strong for the most-bound orbits. As I show below, the least-bound orbits do not strongly lock to the 
black hole's spin until they have very nearly reached the innermost orbit. The NHS orbit's dL z jdi behavior follows 
from the fact that the most-bound orbit locks to the horizon more readily than the least-bound orbit. 

This behavior could have interesting observational consequences. It has been well-understood for some time that 
the inclination angle of an inspiraling body increases due to gravitational- wave emission [^|,^| . Gravitational waves 
carry L z away from the orbit, so that dL z /dt < 0. Since "normal" orbits have dL z /db < 0, it follows that dt/dt > 0. 
If dL z /di switches sign, then dijdi will switch sign as well: an inspiraling body will evolve toward an equatorial orbit. 
If the change in 1 is large, it could have a large effect on the gravitational waveform. For example, the spin-orbit 
modulation of the wave's amplitude and phase is due to motion in the 6 coordinate. This modulation will be reduced 
as the body's motion reduces. 

Since the size of NHS orbits is significant near the "astrophysically maximal" value a — 0.998M (the value at 
which a hole's spin tends to be buffered due to photon capture from thin disk accretion; see Ref . jl0| ) , astrophysical 
black holes might spin quickly enough for NHS orbits to play some important role. This motivates a careful analysis 
to see what role, if any, NHS orbits might play in gravitational-wave sources. I use the code described in Rcf. |J 
(which uses the Teukolsky and Sasaki-Nakamura equations jll],[l2| to compute the flux of energy and L z carried away 
from the orbit by gravitational waves) to study how a small body's motion evolves as it spirals through the NHS 
region. By computing the change (dr/dt,di/dt) at a large number of points, it is straightforward to construct the 
inspiral trajectory for a small body p3| . I find that the total change in inclination angle as a body spirals through 
the NHS region is very small — at most, 81 ~ 1° — 2°. (See Fig. |^; note that the shape of the curves in this figure are 
independent of the values of \i and M, although the timescales strongly depend on fi and M.) 

The inspiral code relies on the fact that circular orbits remain circular as they evolve on an adiabatic timescale 
fj^-|l6||. It is thus explicitly restricted to the evolution of circular orbits, and cannot say anything about the evolution 
of eccentric orbits, which are much more realistic as LISA sources |l^Jl^]. Based on the leading order corrections to 
di/dt seen from a post- Newtonian analysis |l5[ , I speculate that the change in 1 might be much larger when eccentric, 
inclined orbits evolve through the strong field. Verifying this, however, will require a strong-field radiation reaction 
formalism that can evolve generic Kerr orbits: orbits that are inclined and eccentric. Such a formalism may need 
to be based on a local radiation reaction force pj[ , although recent discussion suggests it may be possible to evolve 
generic orbits using gravitational- wave fluxes alone pej] . The possibility that inspiral through the extreme strong field 
may leave an observationally significant imprint on the system's gravitational waveform will hopefully motivate future 
activity and progress on this problem. 

Throughout this paper, the quantities t, r, 9, and <fi refer to the Boyer-Lindquist coordinates. A prime on any 
quantity denotes d/dr. (Note that in Ref. Q I erroneously wrote that prime denotes d/dr.) Section |ll| reviews the 
properties of circular Kerr orbits, providing formulas that are useful for describing their conserved quantities E, L z , 
and Q in the very strong field. These are used in Sec. |H]| to study the N HS orbits. The NHS orbits are developed and 



mapped out as functions of spin and radius in Sec. HlTAL Section III B then re-examines them from the viewpoint of 
the "zero angular momentum observer", or ZAMO. The ZAMO makes local measurements of the orbital properties, 
and can see that their defining behavior dL z /dt > arises due to the locking of NHS orbits to the black hole's spin. 



In Sec. [V, I examine the trajectory of a body that is inspiraling through the NHS region under gravitational- wave 
emission. Here I show the very small change in t as the body spirals in, and argue that eccentricity might impact this 
result greatly. Some concluding discussion is given in Sec. |yj 

II. CIRCULAR ORBITS OF KERR BLACK HOLES 

Geodesic orbits of a Kerr black hole with mass M and spin per unit mass a are governed by the following four 
equations Q: 

£ 2 = [E(r 2 + a 2 ) - aL z ] 2 - A [r 2 + (L z - aE) 2 + Q]=R 1 (2.1a) 
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E 2 (f-) = Q - cot 2 9L 2 Z - a 2 cos 2 0(1 - E 2 ) = 6 2 , (2.1b) 
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The quantities E, L z , and Q ("energy", "z-component of angular momentum", and "Carter constant") specify a 



family of orbits, and are conserved along any orbit of the family. Here, S = r 



'■ cos 2 9 and A = 



2Mr 



Eqs. fl2.1a| ) and (2.1b) have been divided by /x 2 , and Eqs. (2.1c) and (2. Id) by /i (where fj, is the mass of a small 



body in an orbit); E, L z , and Q are thus the specific energy, angular momentum and Carter constant. Also, a > 0; 
prograde and retrograde orbits are distinguished by the orbit's tilt angle rather than the sign of the hole's spin. 

A circular orbit must satisfy R = 0, R' = 0; to be stable, it also must satisfy R" < 0. These conditions are met for 
some set of orbits everywhere outside the innermost stable circular orbit, or ISCO. The ISCO lies at H 



nsco/M = 3 + Z 2 - [(3 - Z x )(3 + Z x + 2Z 2 )\ 
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(2.3) 
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The ISCO varies from r — 6M for a Schwarzschild black hole to r — M for an extreme Kerr hole. 

At all r > risco there exists a family of circular orbits, each member having a different inclination angle t. We are 
interested in parameterizing these orbits as functions of r and i. Consider first the weak-field limit, r 3> nsco- Ryan 
H has provided formulas which, with some manipulation, give E, L z , and Q as functions of r and t: 
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(2.6) 



Note there is a sign error in Eq. (7) of Ref. ||, as can be seen by taking the zero eccentricity limit of Eq. (6) of Ref. 
5fl.] One can then get L z from the definition of the inclination angle: 



cos L 



(2.7) 



Two features of these formulas are particularly noteworthy. First, energy monotonically increases as l increases: 
dE/db > for all parameters. This turns out to be true everywhere, not just in t he weak field . Second, L z 
monotically decreases as l increases provided we don't abuse the applicability of Eqs. ( |2.6|) and (2/7): dL z /di < 
except when r/M < (6a/M) 2 ^ 3 . This range is not even close to the weak-field, so there is no reason to believe that 
this result is at all physically relevant. Nonetheless, it foreshadows the behavior of the nearly horizon skimming orbits. 

Turn now from the weak field to the strong field. As is conventional pp|Jl5|], I will use Eq. ( p.7| ) to define the 
inclination angle^j even in this regime. The most-bound orbitQ (the orbit with the smallest orbital energy) at each 
radius is the prograde, equatorial orbit. Its constants are given by || 
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(2.8) 

(2.9) 
(2.10) 



2 As discussed in Ref. ||, this angle doesn't necessarily accord with intuitive notions of inclination angle. For example, except 
when a = 0, t is not the angle at which most observers would see the small body cross the equatorial plane. 

3 In this paper, the terms "most-bound", "least-bound" and "marginally-bound" describe orbits at each radius. This contrasts 
with the usage in, e.g., Ref. m where these terms refer to properties of all orbits, regardless of radius. 
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where v = y/M/r and q = a/M . At fixed radius, the orbital energy increases as the tilt increases from the most-bound 
orbit at t = 0° to the least-bound orbit. The least-bound orbit's characteristics depend upon r and the black hole's 
spin. If r > r rct , where la] 



r ret /M = 3 + Z 2 + [(3 - Zi)(3 + Z 1 + 2Z 2 )} 



1/2 



(2.11) 



then the least-bound orbit is just the retrograde, equatorial orbit. This orbit has Q — and t = 180°; expressions for 
its energy and angular momentum can be found in Ref. ||. For radii risco < r < r ret, the least-bound orbit is the 
marginally- stable orbit: the orbit which satisfies R = 0, R' = 0, and R" = 0. This orbit has the maximum allowed 
inclination angle i max at that radius. Any orbit tilted at a larger angle is unstable to small perturbations and will 
quickly plunge into the black hole. 

For the rest of this paper, I will focus on the extreme strong field of rapidly rotating black holes. The orbits of 
interest are well inside the radius of the retrograde orbit. Hence, the energy i? lb , angular momentum L l z , and Carter 
constant Q lb of the least-bound orbit will be determined by numerically solving the equations R = 0, R 1 = 0, R" = 0. 

To compute the properties of a circular orbit, pick two of its constants — e.g., the orbit's radius r and angular 
momentum L z — and solve R = = R' to find the other two. This yields the following solution for E(r,L z ) and 
Q(r,L z ) §: 



E(r,L z )= 



Q(r,L,)-- 



2 L 2 z (r- M) + rA 2 



aL z M (r 2 - a 2 ) ± A^r 5 (r - 3M) + a 4 r(r + M) + a 2 r 2 (L 2 z - 2Mr + 2r 2 )} 



[(a 2 + r 2 )E(r,L z )~aL z ]' 



2 E{r, L z ) 2 -2aE{r, L Z )L Z + L 2 Z ] 



(2.12) 



(2.13) 



There is a sign choice in the denominator of Eq. ( 2.12| ). In general, only one choice is physically meaningful at a given 
value of r. The argument of the square root in the denominator of Eq. (2.12) goes to zero at so me ra dius rbranci^a); 
the plus sign corresponds to r > rbranch(fl); and the minus sign to r < ^branchlja)- In Ref. Eq. ( 2.12[ ) was used with 
the plus sign only since the focus in that paper was on comparatively large radius [in all cases, rb ranc h(a) is close to 
2M). In this work, since I will focus on the extreme strong field, both signs are needed. 

Rather than using Eq. (2.12), I will avoid this sign ambiguity by starting with r and E, and re-solving the system 
R = = R! for L z and Q. Then, 



L z (r,E)-- 



Q(r,E)-- 



EM (r 2 - a 2 ) - Ay/r 2 (E 2 - Tf+rM 
a(r- M) 

[(a 2 + r 2 )E-aL z (r,E)] 



[r 2 + a 2 E - 2aEL z (r, E) + L z (r, E) 2 } 



(2.14) 
(2.15) 



(There ex ists a second solution for L z which has a + sign in front of the A, but it is not physically meaningful.) Note 
that Eq. ( |2.14 ) does not behave well as a — > 0. This is because of a degeneracy in this limit: knowledge of any three 
of the parameters E, r, L z , and Q suffices to determine the orbit (because of spherical symmetry). Since this paper 
deals with a ~ M, this issue is irrelevant here. 

Assembling strong-fiel d or bits now reduces to a simple recipe. First, pick an orbital radius 
to vary from £ mb [Eq ( §!§ )] to E lh (found by solving the system R = 0, R' = 0, R" = 0). 
and Q with Eqs. (2.14) and ( 2.15| ). Parameterize each orbit by its inclination angle i [Eq. 
radius. 



Allow the orbital energy 
For each energy, find L z 
7fi]. Repeat at a new 



III. NEARLY HORIZON-SKIMMING ORBITS 



A. Overview 



Consider for a moment the extreme Kerr limit, a — M. From Eq. (2.3), the ISCO is located at r — M, which is 
also the coordinate of the event horizon. It is not difficult to show that there exists a set of orbits at this radius, with 
the parameters 



2M/V3 < L z < V2M , 
E - L Z /2M , 
Q = 3L 2 z /4- M 2 



(3.1) 
(3.2) 
(3.3) 
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These are the horizon skimming orbits. They obey dL z jdi > 0, similar to the behavior seen when the weak-field 
results for L z are pus hed into t he v ery strong field. This is opposite to the typical behavior, as exemplified by the 



correct usage of Eqs. ( |2.6| ) and (|2.7D . 

Focus for now on the most-bound and least-bound orbits. These orbits bound the behavior of all orbits at each 
radius. As described in Sec. D, it is simple to solve for L l z as a function of radius, at least numerically. Doing so for 
a = M, we find that there is a region stretching to r ~ 1.8M in which Ll, b > £™ b . The orbits in this domain have the 
same dynamical characteristics as Wilkins' horizon skimming orbits, so we shall call them "nearly horizon skimming" 
(NHS) orbits. 

Figure |l| illustrates the NHS region, and contrasts it with the "usual" behavior of Kerr orbits. It plots Uf h and 
L l z for black holes with a = 0.95Af (top panel) and a = M (bottom panel). For a = 0.95M, the least-bound and 
most-bound orbits coincide at the ISCO. Moving out in radius, the most-bound orbit's L z grows and the least-bound 
orbit's L z shrinks. This makes sense intuitively, since the inclination angle of the least-bound orbit grows as we move 
away from the ISCO. (Eventually, it tips over completely to i — 180°, and becomes the retrograde equatorial orbit.) 
The lower panel of Fig. [j] shows the behavior when a — M. We see the horizon skimming orbits at r — M and the 
NHS orbits stretching out to r ~ 1.8M. At that point, Is~ and Uf° cross over. Moving further out in radius, they 
behave in the "normal" way. 

There must exist some critical spin value, 0.95M < onhs < Af, at which NHS orbits first come into existence. 
The NHS orbits are bounded by two radii, risco and r cross [where the angular momentum of the least-bound and 
most-bound orbits cross: £™ b (r cross ) = L l z {r cross )]. Decreasing a must decrease the size of the NHS region — nsco 
and r cross approach one another. The spin at which these radii coincide and the NHS region vanishes is onhs- 



Recall from Eq. ( 2.12] ) that there are two solutions for orbital energy, only one of which is usually physical. At 



the crossover point, however, both of these energies must be physical: that with the minus sign (which is larger in 
magnitude) gives E lh , and vice versa. When nsco = Across the two energies likewise coincide, since the most-bound 
orbit is also the least-bound orbit at the ISCO. Thus, onhs is the spin that satisfies 

£+[?*isco(aNHs), £* b ,a N Hs] = E^[r ISCO (a NHS ), Lf° ', a NH s] (3.4) 



[where E+ denotes the two roots given in Eq. (2.12)]. The solution to this is a^ns/M = 0.952412 

Figure || illustrates the change to the NHS region as spin is varied, vanishing altogether when onhs is reached. The 
size is still significant near a — 0.998A/. This is interesting, since a = 0.998A/ is probably the largest spin value that 
we can encounter in nature fllOf. This opens the possibility that NHS orbits may play a role in astrophysical processes. 



B. Nearly horizon skimming orbits as seen by a ZAMO 

The NHS characteristic dL z /di > can be better understood by considering these orbits from the viewpoint of 
the zero angular momentum observer, or ZAMO H], The ZAMO is the observer that corotates with the coordinate 
system, such that its angular velocity as seen at infinity is 

2aMr 

WZAMO = T-J— 2Y2 2 A ' dl ' ( 3 ' 5 ) 

If one imagines spacetime to be dragged into a whirlpool-like flow by the black hole's rotation, then the ZAMO is the 
observer who simply rides along with the flow. (In accordance with this viewpoint, Wald J22| calls the ZAMO the 
"locally non-rotating observer".) 

The ZAMO examines orbits in its local neighborhood. This allows it to interpret the motion of a body in a strong- 
field orbit with special relativistic formulas. For instance, the ZAMO measures the small body to have velocity v and 
energy Ei OC3 i = (1 — v ■ v)^ 1 / 2 =7. This local energy is not a constant of the body's motion. It is related to the 
conserved energy by the formula [[7| 

E = ai^local + WzAMO^z , (3.6) 

where 



a y (r 2 + a 2 ) 2 -a 2 A sin 6» 2 ^ 

is the lapse func tion. Knowing the energy and the angular momentum of the small body then tells us the body's 
speed v = \fv ■ v as seen by the ZAMO. 
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Applying Eq. (3.6) shows that a body's speed is smallest in the most-bound orbit and highest in the least-bound 
orbit, varying smoothly between the two. This is due to the dragging of spacetime by the hole's spin. Let us contrast 
the prograde and retrograde orbits, i = 0° and i = 180°. The prograde orbit is moving "downstream": part of the 
motion needed to keep it in orbit is provided by the dragging of incrtial frames. It can orbit with relatively small 
velocity. The retrograde orbit, by contrast, must "swim upstream": it must overcome the dragging of incrtial frames 
on top of the motion needed to stay in orbit. It therefore is more energetic than the prograde orbit. The angle i 
smoothly varies the orbit between these extremes, so that larger l corresponds to larger energy (and larger speed). 

The velocity of a body in a non-equatorial orbit has components in both the 9 and <j) directions. The 9 component 
goes to zero, however, at the orbit's turning points (# ma x, ^min), when it reverses in 9. (F ormulas for computing 



?max/min can be found in Ref. |9[; they are just the angles at which d9/d,T = [cf. Eq. fl2.1b| )1.) At these two points, 
the velocity is purely along </>, and the body's motion is fully described by the component v i = v ■ e i (where is the 
^-component of the orthonormal basis that the ZAMO uses to make measurements). Evaluating Eq. ( |3.6[ ) at # ma x/min 
and writing i?i OC ai = 7 gives a condition for at the turning points. 

The top panel of Fig. ^ shows ^(^max/min) f° r the most-bound and least-bound orbits at several interesting spins. 
As expected, bodies move quite a bit faster in the least-bound orbit than in the most-bound orbit. Perhaps more 
interestingly, becomes substantially smaller towards the innermost orbits. This is because of the "freezing" of 
physics near the hole's event horizon; see Ref. (particularly Sees. IIC1 and IIIA4) for further discussion. Close to 
the horizon, a body's motion locks to the hole's spin, and it is dragged into rigid corotation. This horizon locking 
causes vi(r) to peak at r ~ 1.5M: the "Keplerian" tendency of a body to move faster as it moves inward dominates 

at large radii (asymptoting to -J M/r in the weak field), but is overwhelmed as the body locks onto the horizon close 
to the hole. 

Consider next the orbit's angular momentum L z . The ZAMO sees L z as the product of a radius of gyration w and 
a locally measured azimuthal momentum p^ ]?J : 

L * = m Pj> = w l v i ■ ( 3 - 8 ) 

The 7 factor causes p^ to be even more strongly peaked than w^; see the lower panel of Fig. [| 

The radius of gyration is a purely geometric quantity. It is just the circumference of the ZAMO's constant r, 
constant 9 orbit, divided by 27r: 

271 



: C(# max / mi n)/27r — — J \J g</></>(r, #max/min) # 
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(r 2 + a 2 ) 2 -a 2 Asin(9 2 , . 

v ' max/mm . „ /„ qn 

• l + a z COS fcr 



T i > ' v ' ' , ~ ' f 

max/ mm 



This function is plotted in Fig. ||. There are no surprises here. In the most-bound orbit, vj increases monotically 
with orbital radius, asymptoting to r at large radius. The least-bound orbit is more interesting: vj is smaller (not 
surprising, since it is a tilted orbit) and is nearly flat as a function of radius, at least over the range of NHS orbits. This 
near flatness is due to the orbit's increasing tilt: #max/min of the least-bound orbit changes such that the circumference 
of the ZAMO's orbit at max /min remains nearly constant. 

Since L z is just the product of the curves shown in Figs. [3] and ^, any unusual features in the behavior of L z must 
arise from features in p^ and vj. Considering these two figures, we immediately see why NHS orbits have L l z > U!f°: 
the least-bound orbit has so much more linear momentum than the most-bound orbit that it compensates for its 
smaller radius of gyration. The linear momentum is so much larger, in turn, because the most-bound orbit is strongly 
locked to the spin of the black hole. The least-bound orbit is also locked for r very close to nsco- However, it does 
not lock as strongly: being so energetic, the least-bound orbit only locks as the very innermost orbits are approached. 
Hence, NHS orbits exist because very close orbits are forced to move in rigid corotation with the event horizon. 



IV. APPLICATION: EVOLUTION UNDER GRAVITATIONAL- WAVE EMISSION 



As discussed in the Introduction, binary systems consisting of a small body spiraling into a massive black hole are 
one of the more anticipated sources of gravitational waves for space-based detectors such as LISA. The NHS region is 
still rather large for a ~ 0.998A/ (cf. Fig. |^), indicating that there might be plenty of time for the properties of NHS 
orbits to influence the gravitational-wave signal of these sources. 
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In the extreme mass ratio, radiation reaction should operate adiabatically. In other words, the timescale for 
gravitational- wave emission to change an orbiting body's parameters (r, E, L z , Q) should be significantly longer than 
an orbital period: T or b/rGw "C 1. (More careful discussion of this point can be found in Ref. |]].) Because the change 
in these constants is very slow, the body's motion is well-approximated as geodesic over small time intervals. It is 
thus useful to regard the body's true, radiatively evolving trajectory as motion through a sequence of geodesic orbits. 
For circular orbits, we regard the body's inspiral trajectory as the evolution of its radius and inclination angle. 

A simple argument shows that a body spiraling through the NHS region should behave rather differently from a 
body spiraling through "normal" orbits. At a given moment, a circular orbit can be represented as a point on the 
(r, L z ) plane. For example, a point on either the top or the bottom panels of Fig. [l] lying between the most-bound 
and least-bound orbit curves represents a physically allowed orbit. Radiation reaction will tend to drive this point 
downward and to the left in this figure: gravitational- wave emission shrinks an object's radius and reduces L z . For a 
"normal" orbit, as in the top panel of Fig. [D, t he body is pushed toward less-bound orbits, increasing l. This behavior 
was predicted by weak-field calculations tfB|]15|] and was recently confirmed in the strong field || . The evolution is in 
the opposite direction inside the NHS region — radiation emission pushes the body toward the sequence of most-bound 
orbits, so that I decreases. 

Decreasing inclination angle could have important consequences for the gravitational-wave signal. Inclined circular 
orbits emit waves characterized by harmonics of = 2tt/T^ and fig — 2n/Tg (where is the period of a full range 
of 4> motion, and Tg is the period of 9 motion) S. Roughly speaking, the 9 motion modulates the gravitational 
waveform. The waves emitted by a body in an equatorial orbit, on the other hand, depend only on fl^ — there is no 
motion in 9, so the amplitudes of fig harmonics vanish. As the orbiting body moves toward the equatorial plane, the 
importance of the fig harmonics decreases, simplifying the waveform. Observing such simplification would be a clear 
signal that the small body is evolving toward the equatorial plane. 

A rigorous analysis of NHS evolution requires a radiation reaction formalism good deep in the strong field. Such 
a formalism, specialized to the adiabatic evolution of circular orbits, is given in Ref. The code described there 
solves the Teukolsky equation |llj to calculate the flux of E and L z to infinity and down the event horizon. It then 
applies the proof that circular orbits remain circular (changing only their radius and inclination angle) when they 
evolve adiabatically under radiation emission to compute the change in Q. From this information, it is simple 

to obtain the change in r and i. By computing dr/dt and di/dt at a large number of points, it is not difficult to 
compute the inspiral trajectory followed by a body that starts in an inclined, circular orbit. An algorithm and code 
for computing such trajectories will be presented at a later date fl3[| . 

Figure [| shows some trajectories in the NHS region for the inspiral of a 1 Mq compact body into a 10 6 Mq Kerr 
black hole. The hole has spin a — 0.998M. The inspiral time for most trajectories is significant. The number of 
orbits around the spin axis of the black hole by the body varies from N<f, ~ 73000 for the inspiral beginning at l = 0° 
to N,/) ~ 400 for the inspiral beginning at i = 40° . The number of gravitational- wave cycles is roughly proportional 
to Nfi, indicating that waves from the NHS region could contribute substantially to a measured signal. The inspiral 
duration scales with M 2 //i, and the number of orbits with M/fj,. These quantities should be of interesting magnitude 
for a reasonably wide range of masses, at least for spin a — 0.998A/. (As the spin of the black hole is dialed down 
to a = cinhSj the NHS orbits lose significance and inspiral through them becomes irrelevant.) Gravitational waves 
produced in the NHS region will form a significant part of the data measured from inspirals into rapidly rotating 
black holes. 

Figure |^ also indicates that the signature of NHS inspiral is not very large. Although t does decrease, as expected, 
the degree of decrease is quite small. For the inspirals shown, the largest change in i is for the trajectory that starts at 
l = 25°: the final inclination is l = 23.8°. The change 8l = 1.2° is a rather paltry effect, not at all a robust signature 
of the strong field. 

It seems that the unique signature of orbital evolution through the NHS region is not likely to be of observational 
interest. Before dismissing the idea entirely, let us speculate for a moment on how eccentricity might affect this 
conclusion. Using a weak-field radiation reaction force, Ryan has calculated the leading order evolution of the 
parameters of an inclined, eccentric orbit of a Kerr black hole due to gravitational- wave emission [j^). His results 
indicate that the eccentricity e can significantly magnify the rate at which l changes: 



This relation is derived under the assumption that r ^> M, so there is no reason to believe it holds in the NHS region. 
However, it could be indicative of the effect that eccentricity has in the strong field. Notice that as eccentricity varies 
over the range < e < 0.8 (a reasonable range for sources that LISA is likely to measure p3fl ), the factor magnifying 
di/dt increases from 1 to about 100. One can speculate that di/dt will be magnified by a similar factor (at least to 
order of magnitude) under strong-field radiation reaction. At this time, there is no way to tell. Hopefully, strong-field 




(4.1) 
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radiation reaction programs 1 19 2(J] will be able to model the evolution of generic Kerr orbits soon. If this speculation 
proves correct, then the gravitational-waves emitted by inspiral through the NHS region will contain a very strong 
signature of the Kerr metric strong field. 



V. CONCLUSION 



The result of this analysis is somewhat equivocal. I have shown that the "horizon skimming" orbits found by 
Wilkins are a subgroup of a larger family of orbits, the "nearly horizon skimming" orbits existing around any black 
hole with spin a > 0.952412M. The signature characteristic of these NHS orbits is that, at fixed orbital radius, 
the z-component of orbital angular momentum increases as the orbit's tilt increases, in opposition to the "normal" 
behavior. This has the consequence that radiation emission, which carries angular momentum away from the orbit, 
tends to drive the system into an equatorial orbit — the system's inclination angle decreases rather than increases. 

The fact that radiation emission drives these orbits toward the equator is a consequence of the motion of bodies 
in the strong field of Kerr black holes. As a body spirals in from large radius, its orbital speed initially increases, in 
accordance with Keplerian intuition that a body in a circular orbit of radius r has a speed v ~ \J M/r. As the body 
gets very close to the black hole, however, its dynamics become dominated by the nearness of the hole's event horizon. 
All physical processes become "locked" to the hole as the horizon is approached. Eventually the motion of orbiting 
bodies synchronizes with the hole's spin. The least-bound orbit does not lock as quickly as the the most-bound orbit: 
close to the horizon, a less-bound orbit moves quite a bit faster than a highly bound orbit, not locking to the hole's 
spin until the very innermost orbits are reached. This effect is shown in Fig. [3| 

As an inspiraling body's inclination decreases, the modulation of its gravitational waveform by the 9 motion 
decreases. Measuring this decreasing modulation could be an observational hallmark of the Kerr strong field. Unfor- 
tunately, at the moment the most mature formalism for computing the effects of gravitational radiation reaction on 
the orbits of small bodies can only tackle the evolution of circular, inclined orbits. In this restricted case, the total 
change in inclination angle as a body evolves through the NHS region is quite small — at most, the orbital inclination 
changes by about 1 — 2 degrees. Although I have not examined the effect of this change on the waveform, it is hard 
to imagine it will be particularly marked. 

It is extremely unlikely that realistic extreme mass ratio inspirals will be circular. When such binaries initially form, 
their eccentricities are likely to be rather close to 1 [|lT| , p~8| , and remain significant by the time that the binary enters 
the sensitivity band of LISA. Weak field calculations suggest that the rate by which the inclinati on a ngle changes 
might be magnified by a rather large factor. If we take the leading order results at face value [cf. Eq. (LI)], the change 
might be so strong that the orbit is driven nearly into the equatorial plane. 

Because the NHS region is so deep in the strong field, however, it is a mistake to infer too much from leading order 
effects. Further progress will require a scheme to calculate radiation reaction in the strong field for generic Kerr orbits 
— orbits that are inclined and eccentric. This may require local radiation reaction forces |l9| , or a clever scheme for 
extracting the change in the Carter constant from the radiation flux p0[ . The fact that observations might be able 
to detect a strong signature of the Kerr strong field will hopefully motivate future progress. 
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FIG. 1. Angular momentum L z for the most-bound and least-bound circular orbits, as functions of radius. The upper panel 
is for a — 0.95M, the lower for a — M. 
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FIG. 2. The region of nearly horizon-skimming orbits for several values of the spin. The solid lines correspond to least-bound 
orbits, dotted lines to most-bound orbits. Notice that the region gets progressively smaller as the spin decreases from a = M, 
vanishing altogether at a — a cr i t ~ 0.952412M. However, the nearly horizon-skimming region remains fairly large at least 
through the vicinity of a = 0.998M, which might apply to some astrophysical black holes. 
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FIG. 3. The azimuthal velocity and momentum at the 6 turning points of an orbiting compact body, as seen by a ZAMO. 
The top panel displays for the most-bound orbit (dotted curve) and least-bound orbit (solid curve) at several spin values; 
the bottom panel likewise displays p^. As expected, a body in the least-bound orbit moves substantially faster than a body in 
the most-bound orbit. Notice that the motion becomes slower as the ISCO (the innermost orbit) is approached. This is because 
the ISCO is close to the horizon at these spin values. Orbits that come close to the horizon become locked to the rotation of 
the black hole. This close to the horizon, locking is substantial. This locking is reponsible for the peaks in these functions: 
moving inward along a sequence of orbits, a body first orbits more quickly, but then slows as its motion locks to the spin. The 
momentum peak is magnified (note the different vertical scales in the two panels) because p^ = 71^ = v^(l — v?)~ 1 '" 2 . 
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FIG. 4. The radius of gyration vj = Circumference(6 l max / min )/27r for the most-bound orbit (dotted curve) and least-bound 
orbit (solid curve). Notice that vj barely varies for the least-bound orbit over this range. This is essentially because the change 
in #max/min compensates for the change in radius such that the circumference of the least-bound orbit remains roughly constant. 
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FIG. 5. The evolution of radius and inclination angle for a small body inspiraling through the NHS region. The small 
body has fi = IMq, the large black hole has M = 10 Mq and spin a — 0.998M. The inclination angle decreases in this 
region, opposite to the usual behavior. The amount of decrease is, however, rather small. [The gap at the end of the inspiral 
trajectories for i < 20° is because of computational limitations: it is computationally expensive to generate a very dense mesh 
of (dr/dt,db/dt) data close to the ISCO.] 
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